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Abstract 

It is well known that a system, S, weakly coupled to a heat bath, B, is de- 
scribed by the canonical ensemble when the composite, S + B, is described by 
the microcanonical ensemble corresponding to a suitable energy shell. This is 
true both for classical distributions on the phase space and for quantum density 
matrices. Here we show that a much stronger statement holds for quantum sys- 
tems. Even if the state of the composite corresponds to a single wave function 
rather than a mixture, the reduced density matrix of the system is canonical, for 
the overwhelming majority of wave functions in the subspace corresponding to 
the energy interval encompassed by the microcanonical ensemble. This clarifies, 
expands and justifies remarks made by Schrodinger in 1952. 

Key words: canonical ensemble in quantum theory; density matrix; typical wave 
function. 

A quantum system in thermal equilibrium at inverse temperature f3 is described by 
the canonical density matrix 

p^ = ^exp(-/?i7(^)) (1) 

where H^^^ is the system Hamiltonian and Z = tr exp{—PH^^^). The usual justification 
for (|T]) is that it is the reduced density matrix of the system when it is weakly coupled 
to a heat bath, and the composite system is described by the microcanonical density 
matrix at a suitable total energy E. 
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More explicitly, one assumes that it is permissible to neglect the relatively small 
interaction between the system and the bath, so that the total Hamiltonian of the 
composite S + B is given by 

H = = + i^(^) . (2) 

The composite + i? is then assumed to be represented by a microcanonical ensemble 
in some energy interval [E,E + 6], where 6 is small on the macroscopic scale, 6 <^ E, 
but large enough for the interval to contain very many eigenvalues. The corresponding 
microcanonical density matrix is 

(3) 

with -P^g g,,.^, the projection to J^e,e+s], the spectral subspace for H associated with 
energies in the interval [E, E + 5] in the Hilbert space = J^^^+^^ = J^'-^^ 
One readily proves (see below) that in the thermodynamic limit, when the size, i.e., the 
number of components N of the heat bath, goes to infinity while E/N — > e, the reduced 
density matrix of the system S 

is equal to pp. Here tr^^^ denotes the partial trace over J^^^^ and (3 = P{e). 

In this note we show how this result can be substantially strengthened: we prove that, 
in the thermodynamic limit, the reduced density matrices of the overwhelming majority 
of the wave functions of S + B are canonical. We call this statement canonical typicality. 
As a consequence of canonical typicality, it follows that the canonical ensemble is even 
more inevitable in quantum mechanics — arising as it does there without the invocation 
of any genuine randomness — than it is classically. Results in this direction were first 
obtained by Schrodinger jH], and later by Gemmer and Mahler [2]; related results have 
been obtained by Tasaki [U]; see later. 

Typicality in quantum mechanics, as well as in classical mechanics, involves a prob- 
ability distribution on the possible microstates of the system, the distribution in terms 
of which "overwhelming majority" is to be understood. In classical mechanics these 
microstates are points in the appropriate phase space and the distribution is then a 
measure on this phase space. To define typicality for quantum systems we shall take the 
microstates to be wave functions, i.e., points on the unit sphere of M' (up to a phase). 
Even with this identification it may not be clear which distribution is appropriate for 
the (composite) system described by the density matrix pe^s. Here we take that to be 
the probability distribution proposed long ago by Schrodinger |31IHj and Bloch [7|: it is 
the (normalized) uniform (surface area) measure ue,5 on the unit sphere in the subspace 
'^[E,E+s\i i-e., the uniform probability distribution over all normalized wave functions 
with energies in [E.,E + 5]. If we expand \1/ in terms of energy eigenfunctions \Ea) 
of if, \E' = 'Y^Ca\Ea), where the sum is restricted to levels with energies in the interval 
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[E, E + 6], then ue,5 corresponds to the uniform distribution on the surface of the sphere 
^|cq,P = 1. This measure was shown in [7] to be, in a well defined sense, the most 
appropriate distribution corresponding to the microcanonical density matrix, i.e., such 
that 

PE,5= [ (5) 



Let denote the reduced density matrix of the system, given that S* + 5 is in a 
pure state G J^[e,e+5], 

= tr(^) |$)(*|. (6) 

We may then ask: for which wave functions \E' is p"^ (approximately) of the canonical 
form We make the standard assumption that both H^^^^ and H^^^ have pure point 
spectrum and are bounded from below. In J^f^^^ (J^^^^) we choose an eigenbasis of H^^^ 
{H''^^), denoted \e[^^) ,\E2^^) , . . . {\e[^^) ,\E2^^) , . . . ), with corresponding eigenvalues 
^ < ^ < . . . (Ep^ < ^ <...). 

First of all we note the following: saying that for the majority of \Ef's, p* is close to 
p^ is equivalent to saying that if is randomly chosen with distribution ue,5 then with 
overwhelming probability, p* is close to pp. From now on we will thus regard \l/ as a 
(Hilbert-space-valued) random variable. 

We begin by recalling the standard derivation of the canonical ensemble from the 
microcanonical. One has for the reduced density matrix p^^^ of the system that 

p(^) = (dim^^,^+,])-i^dim(J^(^))|Ef))(Ef I , (7) 

i 

where 

« Ie-eI^\e-eI'^>+S] ^ ' 

is the spectral subspace for H^^^ associated with energies in the interval \E - Ef^, E - 
Ef^^S\. It is then more or less clear, and can be rigorously proven under suitable condi- 
tions, that when the bath is sufficiently large p^*^-* ~ p/j, with fi = dS{E)/dE where S{E) 
is the bath's entropy. This follows from the basic fact that S{E) ^ logdim(.^]j e+s])^ 
so that 

dim(^(^)) ^ e^(^-^''') ^ e^(^)-^^''' ~ e"^^'"' . (9) 

More precisely, one proves that p^^^ pp in the thermodynamic limit with (3 = ds{e) / de, 
where s{e) = lim[S'(E)/A^] and e = lim[E/iV]. 

Thus to demonstrate canonical typicality it suffices to establish that (|7j) holds, at 
least approximately, when p^^'^ is replaced by p"^ for typical ^ G J^[e,e+5]- A key step 
of our argument is to note that the uniformly distributed random vector can always 
be regarded as arising by normalization 
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from a Gaussian random vector $ G J^[e,e+s] with mean zero and covariance given by 
the identity operator on Jif[E,E+5]- This means that in the decomposition 

t=$:$:c.,iBf'>iBf'>, (11) 

i 3 

the real and imaginary parts, ReCjj and ImCjj , of the coefficients are independent real 
Gaussian random variables with mean zero and variance 1/2 for those i and j for which 
^(S) ^ ^{B) ^^^^-^^ (^^^^ fj,, ^ Q otherwise). We obtain from (HH) that 

^ = Yl\^f^)\^^) (12) 

i 

with 

\<^^) = J2a^Ef^), (13) 
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where Jj is the set of bath levels j such that E^^^ G [E — E^^\ E — E^^^ + 6], whence 
the reduced density matrix (jH)) is of the form 



P 



|$||2 I / \ I 

" «,« 



fS) fS) (S) 

Now, if 5 were so small that the system's energy spacings AE^ = E\j^[ — Ei are all 
greater than 5, then the relevant energy intervals /j for the heat bath would be pairwise 
disjoint and the pairwise orthogonal, 

($,1$,,) = 5.Hi$iir. (15) 

We argue now that equation (fT^ will continue to hold, at least approximately, even 
without the above assumption on AEi . This is so because, when Jj and /j/ have 
significant overlap, the contributions to $j and corresponding to the sum over the 
terms in p3|) belonging to both /j and /j' will typically be approximately orthogonal. 
To see this note that these sums form two independent random vectors, with uniformly 
distributed directions, in a high- dimensional space. As such they are, with probability 
approaching unity, nearly orthogonal since the expected value of the absolute square of 
the dot product of two independent random unit vectors in an n-dimensional space is, 
by symmetry, 1/n. 

From the representation p3|) we have that 

i 

so that is the sum of the Ni = dim(J^^^'*) independent, identically distributed 

random variables |Cjjp with mean 1. It thus follows, by the law of large numbers, that 
typically 

||$,f ^ dim(^('')) . (17) 
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We thus have for the reduced density matrix using (HH), that typically 

p* ^ (dimJ^s,s+,])-i 5^dim(^(^))|Ef , (18) 

i 

which is what we needed to show. 

Concerning ()18|). we remark that it follows merely from the fact that the reduced 
density matrix does not depend upon ^> for typical ^ G ^e,5i that whenever the 
reduced micro canonical density matrix p*^"^) p^, the same is true for p* for typical 
applying the partial trace tr*^-^) to (0), one obtains that p^^^ = J UE,5{d^!)p^ ^ p*, for 
typical ^ . 

Some essential parts of the argument we have presented here have already been 
described by Schrodinger in an appendix, written in 1952, to his book on Statistical 
Thermodynamics We note that Schrodinger in .5] made the assumption that "in a 
big system . . . the amplitude- squares . . . are on the average [in time] equal for ... 
eigenfunctions belonging to the same [energy] eigenvalue . . . ." He uses this assumption 
for the combined system S+B to derive "exactly the same canonical distribution between 
the amplitude-squares, as is in the customary treatment said to indicate the probability 
of the system being on this or that energy level." However, Schrodinger neither connects 
the assumption with typicality, nor his conclusion with the reduced density matrix p*, 
which he does not even mention. His concern is rather with showing that one need not 
regard a system in thermal equilibrium as being in a energy eigenstate. As he states 
in his preface, "To ascribe to every system always one of its sharp energy values is an 
indefensible attitude." 

Tasaki has studied, as we do here, the reduced density matrix of a system coupled 
to a heat bath when S+B is described by a wave function ^> . He shows that for a special 
form of the coupling Hamiltonian the long-time average of p*^*) is canonical. He then 
argues in a heuristic way that also for typical large times, p"^*^*-* ~ pp. (In this argument, 
there is a hidden typicality assumption on the initial wave function \l/(0), namely that 
all of the energy expansion coefficients are small, which is true of most wave functions.) 
His argument does not yield our stronger statement that p**^*^ ^ pp even at t = for 
typical wave functions. Tasaki also includes some examples that he studied rigorously, 
concerning energy eigenstates of the composite or states that are initially product states 
vj/C-s) ^ -^{B) ^-^j^ ^{s) eigenstate of H^^\ 

More recently, Gemmer and Mahler j2] have established canonical typicality under 
the assumption of very large degeneracy for energy eigenstates by computing appropriate 
Hilbert space volumes. 
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Note added: In a very recent preprint, Popescu, Short, and Winter have established 
canonical typicality under great generality by invoking Levy's Lemma 0. 
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